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The issues related to the structure refinement of Pb(Zr a; Tii_ a ;)03 (PZT) solid solutions are dis- 
cussed. Particular attention is paid on the modelling of the co-existing phases in the vicinity of the 
morphotropic phase boundary (MPB), where local symmetry is often significantly lower than the 
average symmetry. According to recent studies, two-phase co-existence in the vicinity of the MPB 
is a thermodynamical necessity. Significantly different structural models for PZT with x in the 
vicinity of the MPB have recently been published. Structural models, based on x-ray, neutron and 
electron diffraction studies, are reviewed and two essentially different approaches were identified: (i) 
a method where space group symmetry was decreased untill the features due to the local distortions 
were 'modelled' and (ii) a method where the highest space group symmetry compatible with the 
powder diffraction data was used together with a model for a local disorder. Related to the method 
(ii), the essential features of a model taking the /ifcZ-dependent line broadening into account are 
summarized. The importance of the cationic disorder is demonstrated by computing the variation 
in Madelung energies due to the variation in structural parameters for PZT using experimental 
data. The underlying theme is to consider how local regions, with symmetry lower than the average 
symmetry, should be modelled. Method (i) often introduces unjustified structural parameters and 
should be avoided. By studying the temperature and composition dependence of the structural 
parameters, including line shape parameters and phase fractions in the case of the two phase sam- 
ples, one can further test the proposed models. The significance of the local structure for electrical 
properties (switchable polarization and electrical conductivity) is discussed by considering the lo- 
cal Pb ion displacements and their dependence on the neighbouring B-cation configuration. The 
connection between the structure and piezoelectric properties, with emphasis on the co-existence of 
the rhombohedral and monoclinic phases, is discussed in light of the recent first-principles compu- 
tational studies. Also problems connected to sample preparation and data collection are pointed 
out. 

PACS numbers: 77.84.Dy 61.12.Ld 61.50.Ah 81.30.Dz 



I. INTRODUCTION 



Lead zirconate titanate [Pb(Zr x Tii_ x )03, (PZT)] solid 
solutions are among the most widely used ferroelectric 
ceramics. The importance of PZT ceramics with mor- 
photropic phase boundary (MPB) compositions is due 
to the exceptionally good piezoelectric properties they 
exhibit 1 , which in turn has motivated numerous studies 
seeking physical mechanisms responsible for this extraor- 
dinary behaviour (for a summary of the various studies 
see ref. 0). Crystal structure studies and the changes in 
atomic positions versus composition, temperature, pres- 
sure or electric field have a central role for the under- 
standing of these materials: nuclei positions, together 
with the electron density, determine the electric polar- 
ization vector on which the numerous applications are 
based on. It is essential to understand how the po- 
larization vector changes as a function of electric field, 
composition, temperature or pressure. Qualitatively, the 
average structural tendencies versus composition, tem- 
perature or pressure of ionic compounds, such as PZT, 
can be understood by considering the Madelung ener- 
gies which decreases with decreasing crystal volume un- 
til the shortrange interactions balance the crystal volume 
to its equilibrium value. For ionic crystals the contribu- 
tion of Madelung energy to the total energy is dominant 



and implies that crystal volume tends to be minimized. 
By tilting the oxygen octahedra the crystal can partially 
compensate for the effect of a larger average i?-cation (Zr 
and Ti) size which increases with increasing Zr content 
x. This is seen from the behaviour of Zr rich samples: 
once crystal adopts Ric symmetry, instead of i?3m sym- 
metry, the volume per primitive cell is smaller. This is 
consistent with the notion that in the case it would cost 
a lot of energy to contract the oxygen octahedra (due to 
the shortrange interactions which cause the mere com- 
pression of octahedra to be energetically unfavourable) 
overall energy gain can still be obtained by tilting the 
oxygen octahedra. This can neatly be expressed through 



the equation Va/Vb ~ 6 cos 2 



1, where Va and Vb 



are the polyhedral volumes for A and B cations and 
(ui) is the mean octahedral tilt angled. This expression 
means that when there is no chance for oxygen octahe- 
dra to decrease its volume the crystal can still be com- 
pressed by decreasing the volume of the cuboctahedra 
surrounding the A cation. Typically the changes in oc- 
tahedral tilts occur consistently versus composition^ or 
temperature^ and thus particular care is necessary be- 
fore anomalous octahedral tilt systems are considered. 
Secondly, in the case of PZT (and other Pb containing 
perovskite ferroelectrics) the crystals suffer from numer- 
ous defects. These defects do have a significant role for 
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ferroelectricity. An example is provided by the Pb ions 
which arc significantly displaced from their ideal sites and 
which respond to the external electric field or pressure 
by adjusting their position with a corresponding changes 
in polarization vectoriii As was discussed in ref. |a, 
it is not clear if Pb displacements are dominantly dy- 
namic or static. Static displacements have commonly 
been modelled by displacing Pb ions towards (110) di- 
rections and using isotropic atomic displacement parame- 
ters (ADP), whereas the dynamic contribution was mod- 
elled by anisotropic ADP parameters. In some cases, 
one can try to distinguish the static and dynamic con- 
tributions by collecting structural data as a function of 
temperatur A In PZT and PbTiC>3 crystals the contribu- 
tion of anharmonic terms to the potential is significant, 
as was demonstrated by Raman scattering studies in the 
case of PbTiOaiSiii and PZTi£, and obscures the sepa- 
ration of static and dynamic contributions. Similar sit- 
uation occurs in many Pb and Bi containing pcrovskites 
and is often believed to be an important factor allowing 
a co-existence of ferromagnetic and ferroelectric order- 
ing in multiferroic perovskitesi^. Related consideration 
holds for Zr and Ti, whose positions are disordered in 
two ways: (i) the probability that the £?-cation site is 
occupied by Zr (Ti) is x (1 — x) and (ii) the fractional 
coordinates of Zr and Ti are different. The confirma- 
tion of the assumption (i) is not straightforward, but 
various powder diffraction and Raman scattering stud- 
ies show that it is a good approximation for most cases. 
By high resolution powder diffraction instruments it is 
possible to study these peculiar features by extracting 
information from the line shapes, in addition to the de- 
termination of the average structure. However, it became 
apparent that a rather complex lineshape is necessary for 
the modelling of the diffraction peak profiles, for exam- 
ples see refs. Q and 0. The structure of the PZT ce- 
ramics with composition in the vicinity of the MPB are 
notoriously tricky to model. In addition to the aforemen- 
tioned structural distortions there are two other reasons 
obscuring the modelling work: (i) in practice, no sin- 
gle phase samples are obtained, instead one observes two 
perovskite phases and (ii) diffraction peak profiles are 
rather complex. Both factors are connected to the fact 
that this 'boundary' separates tetragonal and rhombohe- 
dral phases which do not possess a group-subgroup rela- 
tionship. The experimental observation that monoclinic 
Cm phase separates tetragonal (space group Pimm) and 
rhombohcdral (at room temperature space group is Rim) 
phases clarified the situation: Cm is a common sub- 
group of Pimm and Rim phases&t Correspondingly, 
instead of morphotropic phase boundary one occasion- 
ally refers to morphotropic phase 6,7 . Thus, it is not a 
surprise that the structural disorder is particularly pro- 
nounced in the vicinity of MPB. This is also partially 
due to fact that typically the structure and composition 
of the grain boundaries and bulk are different. This is 
revealed from the Bragg reflection line shapes, which are 
almost 'Gaussian' like at room temperatur e) 7 ! 16 in the 



case of the less abundant phase. Even though it is still 
an open question whether Cm is stable phase in PZT 
system it provides an explanation to the peak value of 
piezoelectric constant via the polarization rotation the- 
ory: first-principles computations demonstrated that the 
polarization rotation provides the lowest free energy path 
along which there is a large coupling between the polar- 
ization and the electric fielcUi, consistently with the ex- 
perimentally observed large electromechanical response. 
Thermodynamical consideration of the polarization rota- 
tion within monoclinic phase and its impact on piezoelec- 
tric properties is given in ref. 0. This study addresses 
the question if it is necessary to invoke the monoclinic 
phase to explain the observed piezoelectric properties. 

Recently, a model which used two monoclinic phases 
(space group symmetries Cc and Cm) was claimed to 
be the correct low-temperature structure of PZT with 
x w 0.53i£, see also note [H and refs. H3 and|M This 
report stated that i?3c symmetry is not a correct choice 
for PZT with x m 0.52 18 . The fact that superlattice 
reflections allow the distinction of R3c and Cc symme- 
tries in favour of the former was overlooked in this study. 
Nevertheless, Cc symmetry predicts reflections well iso- 
lated from those of the i?3c symmetry. To summarize 
the models proposed in refs. [l8| and |2(J, we note that 
the phase which traditionally has been assigned to the 
rhombohcdral symmetry!, was assigned to Cm symme- 
try and space group symmetry PAmm (at room temper- 
ature) or Cc (at low temperature) was assigned to the 
phase which was modelled by Cm symmetry in refs. 5 6 
and [J Thus, significant changes to the well known PZT 
phase diagrami were proposed. To understand the prob- 
lems related to such a model (abbreviated as Cc + Cm) 
and the reason why the use of space group Cc was pre- 
viously rejected^ it is necessary to summarize the crucial 
role of the local disorder resulting in /ifc^-dependent line 
broadeningSS. Erroneous space group symmetry assign- 
ments likely result in once this line broadening is compen- 
sated for by reducing the space group symmetry^. The 
/i/cZ-dependent line broadening is an inherent property of 
PZT powders due to the local distortions, and is already 
seen in tetragonal Ti rich compositions and even at high 
temperature cubic phased. Namely, the Bragg reflection 
widths were anisotropic already at the cubic phase, al- 
though the linewidths were an order of magnitude smaller 
than in the case of the low-temperature phases. It is 
well known that broad peaks are observed in Raman 
spectra collected at high temperature cubic phase, even 
though no first order Raman scattering is allowed for 
the ideal symmetry. For example, see ref. 12a , where 
the modes appearing in cubic phase were assigned to the 
first order scattering, activated by the local disorder, by 
studying the temperature dependence of the intensities. 
The same behaviour was observed in a closely related 
Vb{m x T\^ x )0 3 (PHT) system with 0.10 < x < 0.40, 
where the peak widths of the 00Z reflections were twice 
as large as the widths of hOO reflections^ 6 -. It must be 
emphasized that no signs of symmetry lowering from the 
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P4mm symmetry was observed in this high-resolution 
neutron powder diffraction study. In contrast, light scat- 
tering experiments revealed that both in the case of 
tetragonal PZT and PHT samples deviations from the 
average symmetry were observedSZiSSiS. Namely, the 
number of Raman active modes was larger than the sym- 
metry found by diffraction techniques would allow. Once 
the diffraction patterns were examined it became appar- 
ent that there was no justification to use lower space 
group symmetry. A model based on lower space group 
symmetry would simply result in physically meaningless 
structural parameters. In this case, the essential point 
was that even if one insist to invoke lower space group 
symmetry in order to explain the observed Raman spec- 
tra the changes in bond lengths should be so large that 
also high resolution neutron powder diffraction patterns 
should reveal it. Although one can decrease the aver- 
age symmetry to decrease the residuals, a more realistic 
model is obtained by assuming that (i) the average sym- 
metry over a length scale of a few hundred nanometers is 
P4mm, and (ii) deviations from this symmetry occur in 
a local scale (of a few unit cells) . The essential point here 
is that these deviations are not periodical, at least not 
in a scale probed by x-ray or neutron diffraction. How- 
ever, ion displacements (particularly those of Pb ions) 
are likely correlated in a scale of a few unit cells 2 . It is 
also worth mentioning that usually there are several line 
broadening mechanisms. Especially powders consisted of 
fine crystallites, typically obtained through wet chemi- 
cal methods allowing a preparation of perovskite phase 
at rather moderate temperatures, show large broaden- 
ing due to the small crystallite size, in addition to the 
anisotropic line broadening 30 . Williamson-Hall plot was 
constructed to separate crystallite and anisotropic strain 
broadening effects^ in tetragonal PZT. This study re- 
vealed that the lattice strain is greater along the [001] 
directions than along the [100] directions. On the other 
hand, samples prepared through solid state reaction of- 
ten have an order of magnitude larger crystallite size. In 
the case of PZT with compositions in the vicinity of the 
MPB the line widths are strongly increasing with decreas- 
ing temperature^, demonstrating that the crystallite size 
and shape alone cannot account for the observed line 
broadening behaviour. Similar, although weaker, phe- 
nomenon occurs in the case of tetragonal PZT and PHT 
samples. 

Transmission electron microscopy (TEM) and elec- 
tron diffraction (ED) studies revealed extra reflections 
in rhombohedral PZTs^i. The origin of these reflections 
was more recently studied by TEM and ED together with 
neutron powder diffraction (NPD) techniques 32 . In this 
study superlattice reflections of the type Ri, i?2, Mi and 
Mn (see Table were observed, consistently with ref. 
l3ll Among these reflections only i?2 is consistent with 
R3c symmetry^. We summarize the observations and 
interpretations given in ref. ][32 as follows: (i) M-type 
reflections were accompanied by satellites, (ii) different 
areas in the same grain resulted in changes in the relative 



TABLE I: . Classification of superlattice reflections com- 
monly found in the rhombohedral phases of PZT by TEM/ED 
studies. Tilts refer to Glazer notation—. Table is adapted 
from ref. 13^ . 

Type Reflection Conditions Tilt 



Ri 


\{hkl} p 


h = k = l 


None 


R 2 


H hkl }p 


h^k^l 


a~ 


Mi 


I{0fcZ} p 


k = I 


None 


M 2 


!{ofcz}p 


k ^ I 


a+ 



intensities of the superlattice reflections so that in some 
cases the extra spots completely disappeared, (iii) by col- 
lecting data on samples prepared from single crystals and 
ceramics it was found that the satellites around M points 
are characteristic feature of only ceramics, (iv) by study- 
ing the temperature dependence of the extra reflections 
they could be assigned to the ferroelectric state (no extra 
reflections were observed at cubic phase), (v) octahedral 
tilts and distortions alone are insufficient to explain Ri 
and M superlattice reflections with the intensities ob- 
served (it was demonstrated that very weak superlattice 
reflections could be obtained using unrcalistically large 
oxygen octahedra tilts and distortions), and (vi) NPD ex- 
periments revealed only R 2 reflections, consistently with 
an average i?3c symmetry. To explain these observations 
models based on locally ordered regions presenting an- 
tiparallel cation displacements were proposed (see also 
ref. where structural models based on NPD experi- 
ments are given) and compared with experiments. Finite 
size effect was proposed as an underlying reason for the 
phenomena observed by TEM and ED techniques. This 
is an important point since it suggests that these puzzling 
phenomena, observed by TEM and ED techniques, do not 
occur in the case of bulk ceramics. However, the tendency 
of Pb ions to form four short bonds with oxygen is a fea- 
ture observed through NPD studies^. This tendency is 
also seen from the models constructed to explain the ob- 
served ED patterns 3 ^. It is also worth to note that similar 
conclusions were obtained through Monte Carlo simula- 
tions: M-like superlattice reflections are not present in 
pure bulk crystals but might be locally induced by sur- 
faces over some range of temperature^. 

After a reasonable model for local distortions or lo- 
cally ordered regions is identified a selection of an ap- 
propriate instrument for the (preferably neutron) pow- 
der diffraction experiments to further test the structural 
models is essential. For example, the Cm phase was iden- 
tified by studying the peak split of the pseudo-cubic 110 
reflections 6 through high resolution synchrotron exper- 
iments. This peak split was beyond the resolution of the 
instrument used in ref. and there is no way to recover 
the lost information. Thus, new refinements do not com- 
pensate for the problems connected to the low-resolution 
data collected with insufficient counting time. The role of 
the spatial composition variation (which in practice can- 
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not be eliminated) in the vicinity of the phase bound- 
ary is an old and still a valid explanation for the two 
phase 'co-existences' and also largely explains the /in- 
dependent line broadening. A qualitative model which 
was constructed by studying changes in structures and 
their mutual phase fractions versus temperature and x 
was proposed in refs. 03 and |3 These aspects are re- 
viewed below. In section[H]the variance in Madelung en- 
ergies, corresponding to the variance in experimentally 
determined structural parameters, is given for PZT. The 
idea of this consideration is not to give a qualitative treat- 
ment for the total energy of disordered PZT system but 
to demonstrate the crucial role of disorder by showing 
that rather small variations in structural parameters cor- 
respond to large variation in Madelung energy. Various 
Pb-ion displacement configurations observed in different 
phases of PZTs are addressed in section ITTT1 This section 
reviews the physical origin of Pb-ion displacement and its 
effect on modelling and electrical properties. In section 
IIVI we discuss the role of the spatial composition varia- 
tion on the PZTs with x in the vicinity of the MPB, as it 
is closely connected to the two phase co-existence. The 
main goal of the section IIVI is to provide an explanation 
for the different type of domains by taking into account 
the recent thermodynamical considerations and to show 
that these ideas are consistent with numerous experimen- 
tal observations. Anisotropic line broadening is discussed 
in section|3 Section lVfl reviews different two-phase mod- 
els proposed for PZT in the vicinity of MPB. Section rVlllI 
summarizes the results obtained through bond-valence 
computations and recent computational studies. 

II. MADELUNG ENERGIES FOR PZT SYSTEM 

The Madelung energy was computed by applying the 
Ewald method^i^. The Madelung energy is given by 
the equation Em = J2kk' z kz' k ctkk' / (to r ), where the sum- 
mation runs over all ions fc and k' (with charges Zk 
and , respectively) in a primitive cell (assumed to be 
charge neutral) and eo is the permittivity of free space. 
Charges were assumed to be integer multiples of the elec- 
tric charge e. Following refs. 1351 and l36l the coefficients 
Qfcfc' are written as 

1 \ v 

OLkh' = 7;<V , H(c\xi>k> - Xofc|/r) - 4fc/c/v / 7T+ 
V 

+Tr/(2c 2 s)J2'G{Tr 2 \h{hM)\ 2 r 2 /c 2 ) exp[27rib(/wW) • 

hkl 

(x fc' - x 0fc )r 2 /c 2 ], (1) 

where H(x) = (1/x) erfc(x) (erfc is the complementary 
error function), G{x) = {1/x) exp(x), Skk' is the Kro- 
necker delta function, x/& is the position of the ion k at 
the primitive cell I, h(hkl) is a primitive reciprocal lattice 
vector (indexed by integers h, k, I and defined without the 
prefactor 2tt) and v is the volume of the primitive cell, 
v = sr 3 . Primes above the summation signs indicate 



that in the first (real space) sum the term x;/fc/ = Xofc is 
omitted, and in the case of the second (reciprocal lattice) 
sum the term h(hkl) = is omitted. Number c controls 
the convergence of the two sums, and was selected so 
that both the real and reciprocal lattice sums were of the 
same order of magnitude. The coefficients a k k' do not 
depend on the value of c (dakk'/dc vanishes for all c). 
However, by selecting too small or too large value for c 
results in a situation where either the reciprocal or real 
lattice sum vanishes (within the numerical accuracy of 
the computational routines), leading to erroneous values 
for oikk' ■ We fixed c at it. 

Crystal structure data for PZT ceramics were adapted 
from refs. 5-24 and 37 (high pressure data). To simplify 
the computational task, we assumed that the fractional 
z coordinates for both B cations are the same (though 
they differed in the case of the tetragonal perovskites, 
see Table VI in ref. 0). Thus, for this case the z coor- 
dinate values for the B cations were computed using an 
equation z(B) = x x z(Zr) + (1 — x) X z(Ti) (this ap- 
proximation was used for x = 0.20,0.30 and 0.40 data). 
Similar approximation was also done for the fractional x 
coordinates of the x = 0.52 PZT sample. In addition, we 
assumed that Pb ions were fixed at origin. During the 
Rietveld refinement they were allowed to shift towards 
the (110) directions. Ionic charges were fixed at their 
nominal values (+2 for Pb, +4 for the B cations, and 
—2 for O). An error estimation for the Madelung energy 
was carried out using the standard deviation values of the 
structural parameters. The error estimates of the struc- 
tural parameters were assumed to be independent. It is 
worth to note that the error estimates are specific to the 
structural model. For instance, for the cubic symmetry 
there is only one structural parameter, the a axis value, 
while there are 11 structural parameters for Cm phase. 
This resulted in significantly larger error estimates for 
Cm phase. 

Figure n shows the room temperature Madelung ener- 
gies for PZT ceramics versus x and, for the x — 0.54 sam- 
ple, versus temperature. The Madelung energy increases 
almost linearly with increasing x (when x < 0.50), as is 
seen from Fig. 2] (a) . This is related to the fact that the 
average bond lengths increase with increasing x. How- 
ever, octahedral tilt (i?3c phase) and collective Pb ion 
shifts (Cm phase) correspond to a small energy gain. 
Both are features which were neglected in the present 
treatment of tetragonal PZT. It was interesting to note 
that, in the vicinity of the MPB, Pimm (x = 0.50), 
Cm and i?3c phases (x = 0.52) had almost the same 
Madelung energies: Cm phase had the lowest energy, 
while the energy of i?3c phase was almost the same, and 
approached the energy of Cm phase with increasing x 
(Fig. H(a)). This is consistent with the observed phase 
transition sequence against x. As the difference between 
Madelung energies of Cm and R3c phases in the vicin- 
ity of MPB is small, it is not a surprise that a prepa- 
ration of single phase samples within this composition 
range is very difficult. Also the evolution of Madelung 
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energies versus temperature indicates that the difference 
between the Madeiung energies of Cm and i?3c phases 
is small up to room temperature, see Fig. ^ (b). The 
Madeiung energy of R3c phase at 583 K should be treated 
with caution, since i?3c phase fraction at this tempera- 
ture was very small (correspondingly, the structural data 
were not very accurate) . It is clear from Fig. ^ (b) that 
the Madeiung energy solely does not explain the observed 
temperature dependent phase fractions seen in PZT with 
x = 0.54 sample 2 ^: Cm and R3c phases had almost the 
same Madeiung energies, and it was only at around 583 
K that the difference became significant. Although the 
given treatment is rather crude, it points out that both 
local cation shifts and octahedral tilts are mechanisms to 
be considered. The first one is significant through the 
whole composition range whereas the latter cannot be 
neglected at higher Zr concentrations. 



III. LOCAL DISORDER: PB ION 
DISPLACEMENTS VERSUS x 

The physical origin of the Pb-ion displacement in a lo- 
cal scale is often connected to the 6s 2 lone electron pair 
(L) which makes displaced position energetically more 
favourable (actually one must also consider the 6p states, 
as discussed below) . An insight to the important role of L 
can be obtained by considering the formation of PbTiC>3 
by alloying the litharge phase of PbO and TiC>2 , as shown 
in ref. |23 In the case of Pbi_ a; (TiO) I solid solution a 
TiC>2 group is substitued for a PbOL group and the extra 
oxygen fills the volume of L. Thus, due to L Pb prefers to 
form four short Pb-0 bonds in a pyramidal configuration, 
a tendency which is also seen in different phases of PZT: 
to fulfill this criteria Pb ions are displaced from their av- 
erage positions (although the resulting PbC>4 pyramid is 
no more symmetric). As was discussed in ref. 0, this oc- 
curs in the orthorhombic, tetragonal and rhombohedral 
phases of PZT. In other words this means that Pb ion is 
not located at the centre of the oxygen cuboctahedra, not 
even in the cubic phase. In ref. local structures capable 
of producing the extra reflections of the type ^,0} p 
(Fig. 7) and {§, ±, \} p (Fig. 8) seen in TEM and ED 
studies at approximately the correct intensities are given. 
This picture gives a qualitative idea concerning the Pb 
ion displacements with respect to its nearest neighbours. 
Correspondingly, one must consider the whole oxygen 
and B-cation network in order to determine the over- 
all Pb-displacement pattern. Table |n] summarizes the 
average structures of PZT system with typical structural 
parameters. Fig. [5] shows the commonly assumed Pb- 
ion displacements for tetragonal PZT, PZT with x in the 
vicinity of the MPB (Cm symmetry) and PbZrC>3. Fig. 
121 demonstrates that there are numerous ways to displace 
Pb ions to form short Pb-0 bonds: in PbZrC>3 Pb ions are 
displaced in such a way that an antiferroelectric ordering 
results in in ab plane 39 . Rough but illustrative picture 
for this can be obtained by considering the tolerance fac- 



tor t = (R A +R ) I (V2(R B +Ro)) where R A , Rb and R 
are the ionic radii for the A- and £?-cations and oxygen. 
In the case of PbZrC>3 t < 1 implying that the Pb ions is 
not able to fill the cuboctahedra, whereas Zr fills octahe- 
dra so tightly that it takes a larger share from the total 
volume by tilting the octahedra. In this sense one can 
'construct' PbZrC>3 by (i) considering the large Zr ions, 
which results in a a~a~c° tilt corresponding to distorted 
oxygen octahedra and then (ii) displacing Pb ions inside 
the cuboctahedra in such a way that four short, essen- 
tially covalent, bonds with oxygen (constituting a PbC>4 
'pyramid') are formed. Very similar situation occurs in 
a closely related PbHf03 which is isostructural and has 
the same space group symmetry as PbZrOg^. For clar- 
ity it must be said that the given treatment is meant to 
be illustrative and not to give any causality relationships 
(in ref. the type of octahedral tilting and deforma- 
tion was seen as a way to accommodate the nature of the 
Pb-0 bonding). Now, the Pb04 pyramids alternate in 
the ab plane in such a way that antiferroelectric order- 
ing is formed. Also the Zr ions experience similar shifts, 
although the absolute variation in Zr-0 bonds (at 100 K 
the shortest and largest Zr-0 bonds were 2.41 and 3.105 
A) is not that large as in the case of Pb-0 bonds (at 
100 K the shortest and largest Zr-0 bonds were 2.044 
and 2.205 A). In contrast to Pb ions Zr ions exhibit anit- 
parallel [001] shifts 3 ^. Also the behaviour of the Pb(l) 
and Pb(2) ions (specified in table ITT1 versus temperature 
was found to be different. Namely, at 100 K the Pb(l) 
ion shift was significantly larger that of the Pb(2). Thus 
changes in Pb(l) position were larger than those of Pb(2) 
once the cubic phase was approached. It is interesting to 
note that in the case of Pb(l) three of the Pb-0 bonds in 
Pb04 pyramid were much shorter than the fourth one so 
that Pb(l) is almost threefold coordinated. In this con- 
text it is worth to mention that Pbam and i?3m phases 
do not have group-subgroup relationship nor is common 
subgroup reported to exists between the two phases (in 
contrast to the case of MPB). Interestingly, there is a 
very narrow (with respect to temperature and x) region 
where PbZr03 and PZT with very small amount of Ti 
with Pbam symmetry transforms into ferroelectric rhom- 
bohedral phase at around 500 K. However, the stability 
region of this phase is very sensitive to impurities^ 1 - and 
Pb ion positions are disordered^. Both phase transi- 
tions Pbam — > i?3m and Rim — > Pm3m are first or- 
der transitions (first necessarily by the given symmetry 
arguments) 4 "^. Since the Pb ion displacements from their 
ideal cubic sites are dominant in PbZrOs (so that the Pb 
ion are shifted with respect to the Zr and oxygen ions, 
which remain at their ideal sites) they essentially mea- 
sure the sublattice polarization. Thus it is possible to 
simplify the structural model which in turn allows one to 
write down a simple free energy expansion. 41 for PbZr03. 

In contrast, Pb-ions displacements average out in 
tetragonal PZT so that there is no net polarization in 
ab plane. This is often modelled by assuming that the 
four (xxO) sites, equivalent in the case of PAmm sym- 



TABLE II: Average space group symmetries and representative values for structural parameters reported for PZT system. Note that sometimes the positions of Zr 
and Ti were refined independently. In the case of the rhombohedral phases the structural parameters refer to the hexagonal setting. Customarily R3c symmetry has 
been used to describe the R3m phase by setting e = (correspondingly the c axis is halved). Thus only R3c phase is specified. Lattice vectors and oxygen octahedra 
tilts are given in terms of pseudocubic lattice vectors, xzr refers to Zr content. Deviations from these symmetries are discussed in text. 
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FIG. 1: The Madelung energies (in units of eoA/e 2 ) of PZT crystals per formula unit (a) versus composition x at room 
temperature, (b) versus temperature in the case of the x = 0.54 sample. Em values correspond to P4mm phase for x < 0.50 
and, in the case of the x — 0.52, 0.53 and 0.54 samples, values for Cm (squares) and R3c (triangles) phases are given. The 773 
K value corresponds to Pm3m phase. For x — 0.40 sample, the structural data for 2.7 GPa value was adapted from ref. l37l In 
the case of the two phase samples, the larger error bars correspond to Cm phase, see text. 
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FIG. 2: Pb-ion displacements in (a) tetragonal (space group PAmm), (b) monoclinic (space group Cm) and (c) orthorhombic 
(space group Pbam) phases. Shaded areas indicate the unit cells. 
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metry, are statictically occupied by Pb ion (this is con- 
trasted to the case where Pb ion is located at la site and 
anisotropic ADPs are used) . The structural models used 
in Rietveld refinements take this into account by assum- 
ing that each site is occupied by a quarter of Pb ion. In 
the case of PZT the four sites are not equivalent, but to 
a first approximation, the occupation probability of each 
site does depend on the local S-cation configuration (how 
the eight nearest B cation sites are occupied). Thus, the 
Pb-ion displacements are different in the case of Zr- and 
Ti-rich areas. Cm symmetry allows that spontaneous po- 
larization has a non-zero value in the a-axis direction. In 
the case of ideal symmetry, this corresponds to coherent 
Pb-displacements. However, due to the local distortions 
following spatial composition variation (and changes in 
bond lengths) also other Pb ion displacements (local dis- 
order) must occur. Thus, there is a delicate interplay be- 
tween B-cations, oxygen octahedra tilts (if any) and the 
corresponding Pb-ion displacements. It is not straight- 
forward to say which displacement configuration is en- 
ergetically the most favourable and is most reliably ad- 
dressed through experiments. This is demonstrated by 
PbZrC>3 and PbHfOs: due to the large space available 
for Pb one can construct various hypotetical PbC>4 con- 
figurations (different from the experimental ones, such 
as a one corresponding to a ferroelectric ordering) under 
the constraint that ions have their nominal valences. It 
is not obvious when a change to another structure would 
be favourable. These are aspects to be considered once 
new materials, such as multiferroics, are designed. 

PZT and Ba(Zr x Tii_ :E )03 (BZT) systems provide in- 
teresting case study, since the changes in their electronic 
band structures versus x are well documented. As was 
summarized in ref. l43l there are two main differences 
between Pb perovskites and those containing closed-shell 
A-cations (such as Ba in the BZT system): (i) the band 
gap in PZT system is almost constant, being 3.45 eV 
in PbTiC>3 and 3.72 eV in PbZrC>3 whereas the band 
gaps in BaTiC>3 and BaZrC>3 are 3.0 and 5.0 eV, respec- 
tively and (ii) the presence of shallow Pb 3+ centers. To 
understand these observations, band structure computa- 
tions and electron-paramagnetic-resonance (EPR) spec- 
troscopy studies were carried out for PZT system, re- 
ported in refs. ^3 and El and summarized below. In 
BZT, the valence band edge consists of O 2p states and 
the conduction band consists of pure Ti/Zr d states 45 . 
The band gap increases in BZT with x because the Zr 
4d states lie 2 eV above the Ti 3d states and Ba has 
little effect on the electronic structure because its states 
lie well away from the band gap. In PZT, the valence 
band has essentially the same width and character for 
all x, because this is determined by the Pb s-0 p inter- 
action. Similarly to the case of BZT system, the Ti/Zr 
d states increase rapidly in energy with x. The tight- 
binding computations for cubic PbTiC>3 and PbZrC>3 re- 
vealed that the conduction band minimum is composed 
of the Ti/Zr d states only at low Zr compositions and 
switches to Pb 6p states with increasing Shallow 



Pb 3+ centers corresponds to a local state slightly above 
the valence band maximum. EPR experiments revealed 
that the Pb 3+ center acquires more Pb 6p character with 
increasing x. This was interpreted to be due to the local 
off-center displacement of the Pb ion. Thus, by lowering 
the symmetry of the center the p character is introduced 
into its wavefunction. The Pb 3+ hole trap binding en- 
ergy was found to be between 0.14 and 0.26 eV— . These 
finding were found to be consistent with Raman scat- 
tering studies of PZT, Nd-modified PbTi0 3 (PNT) and 
Nd-modified PZT (PNZT) bulk ceramics using ultravio- 
let (wavelength 363.79 nm which corresponds to 3.4 eV, 
close to the band gap energy) and visible light (wave- 
length 514.532 nm which corresponds to 2.4 eV) 46 . In 
the case of uv light, several high-frequency modes above 
1000 cm -1 were observed. Most interestingly, the fre- 
quency of the mode at around 1170 cm -1 in PZT ce- 
ramics showed anomalous behaviour at the well-known 
phase transitions, corresponding to the changes in Pb ion 
environment. The use of different systems (PZT, PNT 
and PNZT) allowed the assignment of these modes to an 
electronic processes in the Pb 3+ hole traps^i, in contrast 
to Zr and Ti. Once this information is combined with 
the Pb ion displacements observed through NPD stud- 
ies one gains very consistent picture of these systems. 
Now, these observations are important not only for the 
understanding factors affecting on switchable polariza- 
tion in Pb-perovskites, but also provide a physical model 
explaining the activation energies Ea of conductivity in 
PZT. Values between 0.15 and 0.18 eV were calculated 
for the Ea assuming Poole-Frenkel mechanism for the 
conductivity 4 ^. In many ways analogous situation ap- 
pears in Bi-perovskites, such as in (BiSc03) 2; -(PbTi03) a ; 
(BS-PT) alloys—. Namely, the large polarization and 
piezoelectric responses in BS-PT system were found to 
be due to the Bi/Pb 6p-0 2p hybridization 4 ^. 



IV. SPATIAL COMPOSITION VARIATION 

The PZT system is frequently assumed to be a binary 
solid solution of PbTiC>3 and PbZrC>3 without solubility 
gap s, as was summarized in ref. 49. According to ref. 
l49l once a solid solution is formed at high temperature 
(typically above 800°C), the chemical composition can- 
not be changed at low temperatures, but the system can 
have temperature-induced diffusionless structural phase 
transitions. The validity of this assumption has been 
questioned in refs. |5fJ and |5l| and it was reported that 
PZT solutions exhibit a positive enthalpy of mixing that 
suggests a tendency towards immiscibility and phase de- 
composition. Thus, it was proposed^ that miscibility 
gaps replace the MPB and the paraelectric to ferroelec- 
tric transition lines of the diffusionless phase diagram 
given in ref. Q (which would be valid only if the cool- 
ing rate significantly exceeds the diffusion rate of Ti and 
Zr atoms). The analysis was restricted to temperatures 
above 420 Kfil, although very similar conclusions should 
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be valid at lower temperatures. Usually the solubility gap 
increases with decreasing temperature. In the present 
case the solubility gap and finite atomic diffusion rate 
would imply that system minimizes its free energy by 
segregating into two different phases possessing differ- 
ent compositions and crystal symmetries. Although it is 
probable that the equilibrium situation is not necessar- 
ily reached through standard sample preparation tech- 
niques, two-phase co-existence is apparent in the vicinity 
of the MPB. From the microscopic point of view this 
means that attention should be paid on how Zr and Ti 
ions are distributed over the B-cation site. This issue 
was recently addressed in ref. through Monte Carlo 
type model where the mean of the distribution of clus- 
ter volume was analysed as a function of x. This study 
revealed that the sum of the mean cluster size of the 
Zr and Ti clusters is a minimum at x — 0.50, with the 
most likely minimum volume in which an equal number 
of Zr and Ti ions is found being 1 nm 3 (for clarity we 
point out that correspondingly there were many 1 nm 3 
volumes with inequal number of Zr and Ti ions) . Corre- 
spondingly, it was argued that the increase in coherence 
length of the nano-domains, proposed as the mechanism 
for the appearance of the macroscopic monoclinic phase 
close to x = 0.50 in ref. 0, may be a consequence of some 
orderingof the i?-cations. According to the model given 
in ref. |2j, the local structures of the R3m and PAmm 
phases were considered to be monoclinic and Cm phase 
was interpreted to correspond to those composition and 
temperature values for which the local regions had grown 
sufficiently that diffraction techniques see a distinct Cm, 
phase. For the forthcoming discussion, it is important to 
notice that the spatial composition variation is larger the 
smaller the spatial length scale is, as is also apparent from 
Fig. 13 in ref. 0. Our next goal is to estimate the mag- 
nitude of the spatial composition variation in simplest 
terms and to relate it to the experimental observations 
on local and average structure. 

Thus, we assume that the distribution of Zr and Ti 
is random. For discussion purposes, we further assume 
that the spatial composition can be divided into two 
parts: (i) differences between the average compositions of 
grains, consisted of domains and (ii) composition varia- 
tion within a domain. To quantify the meaning of spatial 
composition variation within a domain we assume that 
the distribution of Zr and Ti obeys the binomial den- 
sity function with probabilities p = x and q = 1 — p 
that the B-cation site is occupied by Zr and Ti ion, re- 
spectively (the solubility gap would further enhance the 
spatial composition variation). In principle, ordering ef- 
fects can be taken into account by introducing condi- 
tional probabilities. We believe that the present approx- 
imation provides a microscopic explanation for the two 
phase co-existence observed in the vicinity of the MPB 
and is consistent with the idea that Cm phase serves as a 
transitional phase between rhombohedral and tetragonal 
phases. To see this, we consider the case of PZT with 
x = 0.50: if we divide a large domain (say, with dimen- 



sions of the order of 1 fim) to cubes containing N primi- 
tive cells, roughly 2/3 of the cubes contain 7Vz r = xN Zr 
ions with Np — y/Npq < Nz r < Np + ^/Npq. If we now 
take the cube edge to be 10 nm (typical spot sizes used in 
TEM and ED studies are between 1 and 5 nm), almost 
1/3 of the cubes have x < 0.49 or x > 0.51 (for large 
values of N binomial distribution can well be approxi- 
mated by a Gaussian distribution). Spatial composition 
variation is the probable origin explaining why some do- 
mains are rhombohedral, whereas some are monoclinic. 
In a transition region (compositions between well estab- 
lished rhombohedral and tetragonal regions) Cm phase 
serves as a transition bridge. The idea of an inhomoge- 
neous distribution of Zr and Ti ions is also consisted with 
the TEM and ED observations according to which there 
is a rather large spatial variation in ED patterns as dis- 
cussed above. Surface phenomena are sensitive to local 
distortions. It is the feature of the binomial distribution 
that the variance is largest at around MPB composition 
(j> ~ q) an d is zero only for pure PbTiC-3 and PbZrC>3. 
This is seen from the line widths which increases with 
increasing x when < x < 0.50. This is supported by 
an atomic pair distribution function analysis carried out 
for PZT powders with x = 0.40,0.52 and 0.60^: the 
MPB is a crossover point with maximum disorder. It is 
also important to note that even if the average compo- 
sition of different grains would be exactly the same, the 
spatial variation can correspond to domains with differ- 
ent symmetries. In practice, some composition variation 
in average compositions of different domains and grains 
exists, as the above mentioned solubility gap suggests. 
Also the composition of grain boundaries is expectably 
different from the inner parts. The most important factor 
behind the two-phase 'co-existences' is probably related 
to the differences in average compositions of different do- 
mains: the domains with larger Zr content can entirely 
transform to rhombohedral phase, while it looks that the 
phase transformation cannot be completed in those do- 
mains containing larger amount of Ti but is frozen to Cm 
phase. We note that Cm phase was interpreted to relieve 
the stress which otherwise would be generated due to the 
interacting rhombohedral and tetragonal domains in ref. 
l53l From this point of view it is not a surprise that the 
Bragg reflections of the Cm phase are much wider than 
those of the i?3c phase, and are dependent on temper- 
ature. It is our opinion that PZT with composition in 
the vicinity of the MPB (and similar systems) should be 
modelled as a two-phase system. 

On the other hand, if one could prepare an ideally or- 
dered PZT sample with x = 0.50, Cm phase might not be 
stable. To see this, we note that in the case of ordered 
PZT with x — 0.50 the Zr/Ti distribution is homoge- 
neous in the primitive cell scale and any deviation from 
this order results in less homogeneous distribution. Cor- 
respondingly, there are no Zr rich or deficient areas and 
in this sense also the mechanism favouring Cm would be 
missing: if there are no local distortions due to the spatial 
composition variation, there likely are no domains with 
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different symmetries (at least not so abundant as the ex- 
periments reveal), which in turn implies that there is no 
need for a phase serving as a bridge between rhombohe- 
dral and tetragonal phases. This might be important for 
a very small grain size powders where most of the grains 
would be consisted of one domain. Experimentally the 
two extrema cases, perfect disorder and order, can eas- 
ily be distinguished by x-ray or neutron diffraction tech- 
niques. As the studies dedicated to double perovskites 
show, ordering is clearly seen from the superlattice reflec- 
tions in x-ray (for instance, see ref. |^] where polymerized 
complex method was used to prepare Ba2MnWOg dou- 
ble perovskites) and neutron powder diffraction patterns 
(see, e.g., ref. HQ). Also Raman spectra are consisted 
with the observed diffraction patterns^S^i. As the nu- 
merous studies show the difference in ionic radii (and, to 
a lesser extent, charge) does not provide sufficient driving 
force for a formation of the double perovskite structure. 
As proposed in ref. 0, partial ordering may occur, al- 
though it is not easy to determine the degree of ordering 
in the present case. Even though it may not be possible to 
prepare double perovskite PZT samples the problem can 
be approached by computational techniques. Finite tem- 
perature Monte Carlo simulations carried out for PZT 
with a broad range of x and random distribution of Zr 
and Ti in the S-cation site yielded a phase transition 
sequence Pimm — > Cm — > i?3m£2 (no interaction be- 
tween different domains was incorporated in this study). 
Virtual crystal approximation (VCA, where Zr and Ti 
ions are replaced by the same fictitious average atom) 
does not reveal a Cm phased. It is an open question 
if VCA could reproduce Cm phase when the interaction 
between domains is taken into account. After all, Cm 
phase always co-exists with rhombohedral phase. The 
two-phase co-existence in the vicinity of the MPB has 
proven to be very puzzling and sometimes one can find 
reports according to which an 'ideal' sample preparation 
route could yield single phase samples. Recently, an equi- 
librium phase diagram satifying the Gibbs phase rule was 
given in ref. OH} The essential result was that, by the 
Gibbs phase rule, the single-phase fields on the phase di- 
agram must be separated by two-phase regions (and Cm 
phase is by no means an exceptional. According to ref. 
\VA the flat energy surface near the rhombohedral phase 
may be the most important factor for high piezoelectric 
constant materials. As the co-existence suggests, this 
seems to be the case in PZT system. 



V. ANISOTROPIC LINE BROADENING. 

The anisotropic (/ifcZ-dependent) line broadening was 
ascribed to a 'microstrain' in ref. [24|. There are several 
factors which together are responsible for the anisotropic 
line broadening: (i) Zr substitution for Ti creates local 
strains due to the size difference, (ii) spatial composi- 
tion variation: Zr-rich clusters take larger volume per 
formula unit than Ti-rich clusters, and (iii) Pb-ion shifts 



depend on the i?-cation environment (i.e., Pb-shifts vary 
spatially as a response to the spatial variation of the 
-B-cations). Although the given mechanisms are inter- 
weaved, the given division serves as a way to estimate 
which one is dominating at different composition regimes. 
Thus, the changes in Bragg reflection intensities and dis- 
placement parameters in the case of small Zr/Ti sub- 
stitutions in PbTiOa/PbZrOs can largely be understood 
by considering the mechanism (i). Thus, following the 
treatment given in ref. l59t small Zr substitution for Ti 
in PbTi03 might be treated by considering the concen- 
tration of Zr atoms x, the difference between matrix (Ti) 
and solute (Zr) atoms AR. Now, the degree of distortion 
is characterized by the mean square atomic displacement 
(U 2 ) — jx(AR 2 ), where j is a constant. Static atomic 
displacements cause a reduction in diffraction peak inten- 
sity and increase in diffuse scattering and can be treated 
in the same way as thermal vibrations. However, this 
treatment is no more adequate for intermediate values of 
x since there is no way of selecting such a matrix that 
experimentally observed features could be explained. Al- 
though the mechanism (i) can explain changes in inten- 
sities (matrix and solute atoms have different scattering 
lengths and cross sections) and diffuse scattering (seen 
as an increased intensity in the Bragg reflection tail re- 
gions: the peak width at half maximum is practically 
not affected by this mechanism) for small values of a;, it 
does not provide a model for changes in Bragg reflection 
widths. This is why the phenomenological model given 
in ref. [lj was adopted: it also takes into account the 
/ifcZ-dependent Bragg reflections widths by introducing 
a microstrain broadening Y 2 S = X^hkx SnKLh H k K l L , 
H + K + L = 4. Laue symmetry imposes restrictions 
on the allowed Shkl terms. Expressions for T| for each 
Laue symmetry are given in ref. l60t 

In addition, we note that rather strong intensity is 
commonly observed between the hOO and 00Z with h — I 
Bragg reflections, which was assigned to locally disor- 
dered regions in the vicinity of the domains walls in ref. 
7. Especially at room temperature it is practically impos- 
sible to assign a precise structural model for this phase, 
since it is not well crystallized. It might also be a pre- 
cursor phase which at higher Zr contents evolves into 
a rhombohedral phase. This interpretation was given 
to the data collected on PHT powders, see Fig. 4 in 
ref. |2(J Nevertheless, scattering from this phase adds 
intensity to certain diffraction peaks, which in practise 
also contributes to an asymmetric line broadening (as it 
is not possible to unambiguously separate it: it is the 
sum from different contributions what one observes in 
the case of powder diffraction). For example, the larger 
and smaller d-spacings sides of the pseudo-cubic 200 and 
002 reflections, respectively, seem to gain intensity. Par- 
ticular care is necessary once this phase is modelled and 
generally one might do better by collecting data on sam- 
ples with different Zr contents and at different tempera- 
tures. 

Although a reasonable structure refinements for 
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FIG. 3: Left panel: Schematic illustration of the way how the 
neglect of anisotropic line broadening affects the position of 
the rhombohedral peaks: once the rightmost peak broadens, 
it pushes the rhombohedral peak (middle) towards smaller 
d-spacing. Righ panel: Observed and computed high resolu- 
tion neutron powder diffraction profile collected on PZT with 
x — 0.53 at 4 K. Pseudo-cubic 200 reflection region is shown. 
Note that the 002m peak (Cm phase) at around 2.06 A is sig- 
nificantly broader than the 220m (Cm phase) peak at around 
2.02 A. Now, the position of the rhombohedral 204_r peak de- 
pends on the way the line broadening is taken into account. 
Figure adapted from ref. [(H] 



medium resolution data collected on PZT samples with 
x < 0.50 were obtained using a lineshape which ignored 
anisotropic line broadening, the situation was quite dif- 
ferent for high resolution facilities, particularly once the 
data were collected on two phase samples. In such a case 
it was essential to use an appropriate profile function. 
Now, when the anisotropic line broadening is neglected, 
the fit in the case of certain weak peaks becomes slightly 
worse (as was observed to be the case of weak superlattice 
reflections, which are less weighted in the refinements). 
The reason for this is illustrated in Fig. [J] in order 
to improve the fit corresponding to the strong pseudo- 
cubic 200 reflections, during the refinement the posi- 
tion of i?3c reflection is shifted toward higher <i-spacings, 
which in turn resulted in a small shift of the weak peaks 
at 2.44 and 1.06 A. The anisotropic line broadening was 
not limited to the d-spacing at around 2 A, and simi- 
lar mechanism was seen at other d-spacings. A struc- 
tural model which reduces the symmetries to 'model' the 
/ifcZ-dependent line broadening takes local disorder into 
account incorrectly, since in this way the disorder is as- 
signed to be periodical, obeying the space group symme- 
try. 

VI. TWO-PHASE MODELS 

A. Pb(Zr .52Tio.48)0 3 and Pb(Zr .53Tio.47)0 3 
samples 

In the context of low-temperature phases the role of 
spatial composition variation was discussed in refs. and 
|24L The room- and low-temperature crystal symmetry of 
Zr rich PZT ceramics is i?3c 4 24,62 , except for the compo- 
sitions with x m 1—. The room temperature symmetry 
of PbZrC>3 is Pbam*2£^. Since spatial composition vari- 



FIG. 4: The consequences of the spatial composition varia- 
tion (at fixed temperature and pressure) in the case of two 
samples with average compositions x\ and X2 in the vicinity 
of (a) the phase boundary separating R3c and Cm phases (at 
X\) and (b) two phase boundaries separating R3c and Cc (at 
X2) and Cc and Cm phases (at X3). Figure adapted from 
ref. EH 



ation cannot be completely eliminated, there must exist 
two phases in the vicinity of MPB. Fig. Q|(a) illustrates 
the consequences of spatial composition variation in the 
vicinity of MPB at low temperature. To allow the ex- 
istence of Cc phase necessitates that there should be a 
narrow region in the x — T plane were this phase is stable 
or metastable. Now, if we were to explain the existence of 
Cc phase, two phase boundaries located somewhere be- 
tween 0.52 < x < 0.54 should be assumed to exists, see 
Fig. |H(b). This in turn, once the spatial composition 
variation is taken into account, leads to three phase 'co- 
existence'. The simplest Cm+R3c model (corresponding 
to Fig. El (a)) was preferred in refs. and l24l since it was 
able to explain the experimental observations in simplest 
terms. At low temperature (4 K24 or 10 K£) the phase 
fraction of the Cm phase was monotonically decreasing 
with x increasing from 0.52 to 0.54, which in turn implies 
that two-phase 'co-existence' is predominantly due to the 
spatial composition variation^. 

It is worth to point out that the model proposed in ref. 
[13 was ruled out also in a recent paper, see the last para- 
graph in ref. [6|J After all, the only difference between 
Crri + R3c model and Cm + Cc model proposed in ref. |6^ 
is that R3c symmetry was used (corresponding to two lat- 
tice parameters and four atomic coordinates) instead of 
the Cc symmetry (corresponding to four lattice parame- 
ters and twelve atomic coordinates, although constraints 
were used to decrease refinable parameters in ref. l6o) . 
Now, Cc is a subgroup of Ric space group. Since no su- 
perlattice reflections characteristic only to Cc phase were 
observed (such should appear at large (i-spacing region, 
if symmetry is lowered from i?3c to Cc), R3c symme- 
try was favoured in ref. |5j. This allowed considerably 
simpler structural model: one is tempted to ask how 
much can be gained by first decreasing the symmetry 
from i?3c to Cc and then decreasing the number of refin- 
able parameters by introducing constraints. This issue 
is discussed on a more general ground in ref. l67l In ref. 
log the use of Cc phase was largely based on the TEM 
and ED observations^ according to which superlattice 
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reflections, which were inconsistent with i?3c symmetry, 
were observed at low temperature. This observation is 
consistent with the observations reported in ref. |22J but 
the interpretation is entirely different. Now, it would be 
interesting to consider: (i) similar possibilities as were 
given in ref. (summarized above), and (ii) if the ob- 
served features are only characteristic to TEM samples 
(layers with thickness between 5 and 100 nm were stud- 
ied and it was further stated that the size distribution 
of Cc phase ranged between 3 to 10 nmS. This is con- 
trasted by the observation according to which the line 
widths of i?3c phase were much narrower than those of 
Cm phased. Neither of the phases was in a form of 
'nanoparticles'. It would also be interesting to see how 
large octahedral tilts are necessary to 'reproduce' the ED 
patterns, including those superlattice reflections which 
were taken as an evidence for Cc symmetry, since they 
were not observed by NPD studies. It is an often over- 
looked issue that Cc + Cm and Cm + Cc models, given 
in rcfs. [H and respectively, are completely different 
(and mutually exclusive): they present models where Cm 
and Cc symmetries were swapped. Although this swap- 
ping might not be so obvious from studies solely based on 
TEM/ED techniques it is clearly seen from X-ray pow- 
der diffraction and NPD models. In the case of TEM/ED 
techniques very small volumes are studied in a time (this 
can also give wrong idea about the phase purity, i.e., if the 
sample with macroscopic dimensions contains one or two 
phases) and often no attempts to determine sufficiently 
accurate lattice parameter were done. Whatever is the 
reason for overlooking this issue, this type of space group 
swapping results in entirely different structural parame- 
ters and phase fractions. Needless to say, this has nothing 
to do with composition variation or differences in crystal- 
lite sizes, see also related discussion in ref. |6^ Despite 
the increased number of refined parameters, the resid- 
uals were still rather high and the differences between 
different models compared in ref. were marginal. For 
instance, R exp was lowest for the Cm + i?3c model which 
was rejected in ref. 0. As discussed in the next section, 
the model given in ref. does not assign octahedral tilts 
to a correct phase. In constrast, the models proposed in 
refs. IBI andl6a assigned octahedral tilts to the same phase. 
In the case of the data shown in ref. the intensity of 
the peak(s) at around 1.06 A was barely above the noise 
level and the Cc + Cm model assigns more reflections 
than there are data points in this region. Thus these 
reflections do not provide justification for the use of Cc 
phase. 

In addition, some reports claimed that dielectric mea- 
surements gave support to Cc phase (for example, see 
ref. I20I) . Although these types of measurements might 
be useful once phase transitions in single phase samples 
are studied, not too much emphasis should be put on data 
collected on multiphase samples. In this context we note 
that the samples studied in refs. and |n] contained 
a nonidentified impurity phase (s), as was revealed by 
the peaks at around 28 and 35 two-theta degrees (corre- 



sponding to d spacings 3.18 A and 2.56 A, respectively), 
in addition to the aforementioned two-phase co-existence. 
This implies that the composition was not well known 
and thus the possibility of compositional and structural 
inhomogeneities should be kept in mind. Although this 
non-perovskite phase was well resolved, it was not in- 
cluded in the Rietveld refinement model considered in 
ref. l70l which in turn results in an error in the structural 
parameters of the perovskite phase(s). This is likely re- 
lated to the fact that the diffraction pattern shown in 
ref. [?3 (reported to have x — 0.52) is reminiscent to the 
diffraction pattern with x — 0.53 shown in ref. 5. In con- 
trast, the diffraction pattern with x = 0.52 given in ref. 
IT^I and the diffraction pattern of the x — 0.52 sample 
reported in refs. |6^ and 00 were reminiscent. 



B. Pb(Zro.54Tio.46)0 3 sample 

The anisotropic line broadening and the absence of the 
superlattice reflection evidencing Cc symmetry in PZT 
powders with x = 0.52 and x = 0.53 provided a mo- 
tivation to carry out a subsequent study using a high 
resolution NPD instrument. To model the peak pro- 
files, GSASSi lineshape 4 by Stephens^ was used in this 
study. In addition, Zr content was slightly increased to 
x = 0.5424 (although it was possible to fit all reflections 
using the Cm+R'Sc model also in the case of x — 0.52 and 
x = 0.53 samples). This allowed a more reliable refine- 
ment and symmetry identification to be done by studying 
the changes in phase fractions versus temperature. Im- 
portantly, for these compositions it has been found that 
oxygen octahedra tilts increase with increasing x (see ref. 
0) and decreasing temperature (this feature is discussed 
in refs. and0). We also note that previously Cc was 
proposed to be a space group symmetry corresponding 
to high isotropic pressure^, whereas the present review 
concentrates on the determination of space group sym- 
metries versus composition and temperature at ambient 
pressure. Now, PZT sample with x = 0.54 provided a 
test for clarifying which phase is the preferred one at 
low temperature. This is of particular interest also from 
the point of view of ah initio computations dedicated 
for PZT according to which the largest piezoelectric ^33 
coefficients are found to be large namely in the rhom- 
bohedral side of the MPB: much smaller values were 
found in the tetragonal side of the MPBZ2. These ob- 
servations are consistent with the experimental observa- 
tions on PZTS and Pb(Zn 1/3 Nb 2/ 3)03-PbTi03 (PZN- 
PT) and Pb(Mg 1/3 Nb 2/3 )03-PbTi03 (PMN-PT)24 (the 
two latter systems also possess MPB separating rhom- 
bohedral and tetragonal phases). The NPD data in 
ref. |2J suggests the following phase transition sequence 
with decreasing temperature in the vicinity of MPB: 
PmSm — > PAmm — > Cm — > Rim — > i?3c where adjacent 
groups have a group-subgroup relationship. Due to the 
composition variation some domains undergo the same 
phase transition sequence at higher (Zr richer) or lower 
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(Ti richer) temperatures. This sequence demostrates the 
role of Cm phase linking PAmm and i?3m phases. The 
co-existence of the rhombohedral and Cm phase is con- 
sistent with the idea that they are energetically almost 
as favourable and rather small changes in temperature, 
pressure or an applied electric field can transform Cm 
phase to rhombohedral phase. First-principles compu- 
tational study revealed that the electric field induced 
phase transformation sequence with an increasing ap- 
plied electric field was found to be reversed when the 
field was released 7 ^, whereas the same was not found ex- 
perimentally in ref. l75i This different behaviour was 
attributed to defects^: the model used in ab initio com- 
putations was defect free, whereas in practice they cannot 
be completely eliminated. In this context it is interest- 
ing to note the domain switching in Pb(Zro.4gTio. 51)03 
sample, where tetragonal and rhombohedral (volumet- 
ric phase fractions 79 and 21 %, respectively) phases 
were present, studied through in situ uniaxial compres- 
sion experiments^. Data were analysed via Rietveld re- 
finement which allowed texture and lattice strains to be 
simultaneously determined. It was found that the rhom- 
bohedral phase responds more readily both to electric 
fields and applied stress. 

As is seen from Fig. 1 in ref. |24| once the crystals 
are compressed with decreasing temperature, the i?3c 
phase was favoured and Cm phase did not transform 
to Cc phase, in strong contrast with the model pro- 
posed in ref. Il8f) . We note that the volume of the i?3c 
phase grows with decreasing temperature (for instance, 
see Fig. 1 in ref. [241). Thus, no evidence was found 
for the PAmm — > Cm — > Cc phase transition sequence. 
When the oxygen octahedra had almost no chance to 
contract, they were rotated to opposite directions along 
the pseudo-cubic 111 axis (tilt system a~a~a~). This 
tilt results in a decrease of the volume of cuboctahe- 
dra around Pb ions. Changes evidencing Cc symmetry 
(and the corresponding tilt system) were not observed. 
Indeed, the phase fraction of Cm phase (which was as- 
signed to Cc symmetry in ref. \ld(l decreased with de- 
creasing temperature, whereas i?3c phase fraction signif- 
icantly increased with decreasing temperature. This can 
be confirmed by a naked eye observation of the diffrac- 
tion patterns in Fig. 1 in ref. Wa\ although the results 
of refinements were also given. The intensity of the peak 
at around 1.06 A was increasing with decreasing tem- 
perature, which further confirmed that its origin is the 
phase assigned to i?3c phase. Also the peak at around 
1.06 A was well fit by i?3c symmetry. In other words, in 
ref. Cc phase was used to model the Bragg reflection 
from the pseudo-tetragonal phase and also the superlat- 
tice reflections. Thus, Cc symmetry was used to model 
Bragg peak from two distinct phases, which is not correct. 
However, the temperature and composition dependence 
clearly show that when the phase fraction of Cm phase 
decreased, the intensity of the superlattice reflection was 
increasing together with the phase fraction of i?3c phase. 
This behaviour is consistent with the composition vari- 



ation: the distortion (monoclinic or rhombohedral) and 
the volume of the distortion depend on average compo- 
sition and temperature. 



VII. CORRELATION BETWEEN THE c/a AXIS 
RATIO AND PIEZOELECTRIC CONSTANTS 

It is of interest to compare the recent computational 
results with the structural parameters obtained through 
structural studies. According to the computations car- 
ried out for PbTi0 3 and Pb(Zr . 5 oTio.5o)0 3 (where Zr 
and Ti atoms were alternatively occupying the £?-cation 
layers perpendicular to the c axis) the fixed volume 
piezoelectric constants 633^ = 633 — 0.5(e3i + 632) and 
ei3,D = ei3 — 0.5(en + era) become very large in the 
range 1.0 < c/a < 1.02 77 . Experiments have shown that 
the c axis remains almost constant for < x < 0.50, 
whereas a axis increases almost linearly for these compo- 
sitions, corresponding to a steady decrease of the c/a axis 
ratio^. However, in the MPB region the c/a axis ratio de- 
creases dramatically: averaging monoclinic a and b axes 
the c/a is found to be 1.01 for x — 0.542i. Fig. shows 
the room temperature c/a axis ratios. It is seen that 
the ratio strongly decreases once x increases from 0.50 to 
0.54. This is consistent with the idea that Cm phase con- 
tinuously transforms to rhombohedral phase. For pure 
PbTiOs ab initio computations show that PAmm phase 
is energetically more favourable than the R3m phase, 
consistently with the experimental observation. This 
means that the electric field, which in practice can be 
applied (without electric breakdown), does not induce 
a phase transition sequence PAmm — > Cm — > Rim 
and, via the polarization rotation model, piezoresponsc 
remains much lower than in the case of PZT in the 
vicinity of the MPB. However, recent computational 
study demonstrated that pressure can induce a MPB in 
PbTiO^. It was found that pressure induced a phase 
transition sequence PAmm —f Cm — + i?3m — * Pm3m, 
together with large piezoelectric coupling coefficients in 
the transition regions 7 ^. 



VIII. BOND- VALENCE SUMS AND ION 
VALENCES 

Computations of bond valence sums (BVS) provides a 
way to test whether the bond lengths are reasonable. Al- 
though it is not sensitive for minor changes in symmetry, 
unless bond lengths are altered, it has proven to be useful 
for confirming the effect of the proposed Pb, Zr and Ti 
ion displacements^. Namely, if Pb ions are constrained 
to be located in their ideal positions in the case of PAmm 
phase, the Pb ion valence would be « +1.8, whereas it is 
quite close to the nominal valence +2 once Pb ions are 
allowed to be displaced towards (110) directions. This is 
actually rather general tendency for Pb ion and occurs 
in rhombohedral phases. This is quite plausible since the 
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FIG. 5: Experimentally determined room temperature c/a 
axis values for PZT ceramics. Lattice parameters were 
adapted from ref. IH (PbTiOs) and refs. and|24| (PZT). 



cuboctahedra is so large that if Pb ion would occupy its 
centre, it would result in valence deficient Pb ions. Thus, 
by forming four short bonds with oxygen more reason- 
able valence values for Pb are achieved. Similarly, the 
average iJ-cation valence v(B av ) = aru(Zr) + (1 — x)v(Ti) 
was +4 when Zr and Ti ions were allowed to have differ- 
ent fractional coordinates. In the case of the tetragonal 
PZT it was demonstrated that if one insists to constrain 
the fractional coordinates to be the same, no Wyckoff lb 
position in the unit cell would correspond to the nom- 
inal valence +4. It also worth to point out that if Pb 
ions are constrained to be at the la site symmetry po- 
sition Pb ADPs are very large. On the other hand, by 
constraining the B cations to have the same fractional 
coordinates yielded negative ADP values. Allowing Pb 
ions to be displaced towards (110) directions and Zr and 
Ti ions to have different fractional z-coordinates resulted 
in physically reasonable ADP values. From this point of 
view it was important to check the validity of the struc- 
tural models by confirming that both anomalous valence 
and ADP values can be eliminated by the same struc- 
tural model. Oxygen valences in PZT systems were close 
to their nominal valence value —2. Also oxygen ADP 
values were reasonable. Thus, as far as NPD studies are 
considered, it is a decent approximation to assume oxy- 
gen octahedra to fullfil the requirements of space group 
symmetry, in contrast to the case of cation positions. 
This also implied that it was not necessary to invoke the 
plausible mechanism, the expansion of Zr octahedra and 
contraction of Ti octahedra, to correct the anomalously 
large Zr and small Ti valences in the structural models 
used to model the NPD data. 

BVS were found to be consistent with the observed 
phase transition sequences so that reasonable valence val- 
ues for each ion were obtained also in the case of the Cm 
and i?3c phases. The relationship between the oxygen 
octahedra and cuboctahedra volumes and octahedra tilts 



were pointed out in ref. .3* Namely, i?3c symmetry al- 
lows oxygen octahedra tilt, which means that the volume 
of oxygen octahedra increases with increasing tilt angle. 
This in turn allows the valence of each ion to be close 
to their nominal values. The same is not possible in the 
case of R3m symmetry and thus the observed oxygen oc- 
tahedra tilts versus temperature and x seen in the case 
of PZT samples with x = 0.52 and x = 0.53 are consis- 
tent with the 'constraints' set by nominal valences^. In 
many ways similar findings, based on density functional 
theory computations carried on PZT modelled with large 
supercells, were reported in ref. |8(J. Perhaps the most es- 
sential difference between the typical Rietvcld refinement 
models and supercell models for PZTs is that the space 
group symmetries used in the latter approach allow dif- 
ferent size oxygen octahedra for Zr and Ti, whereas they 
are constrained to be the same in Rietveld refinements 
(one actually replaces both Zr and Ti by a 'pseudo-atom' 
whose scattering length (NPD) or cross section (X-ray 
diffraction) is a weighted average over Zr and Ti ions). 
This explains why the individual bond-valences for Zr 
and Ti were found to be close to their nominal values in 
the case of the supercell computations^. In this sense it 
is surprising to see that both approaches support the idea 
that oxygen position do not suffer from local distortions 
to the extent cation positions do. 



IX. CONCLUSIONS 

The pitfalls related to the modelling of the crystal 
structures of PZT ceramics with different compositions 
and at different temperatures were reviewed. The role of 
the spatial composition variation, resulting in two-phase 
'co-existence' in the vicinity of the phase boundaries and 
anisotropic line broadening, were reviewed. The condi- 
tions set to sample preparation, powder diffraction in- 
strument and line shape were addressed. Focus was put 
on the structure analyses on PZT ceramics with composi- 
tions in the vicinity of the morphotropic phase boundary 
(MPB). Recent studies pointed out that the two-phase 
co-existence in the vicinity of the phase boundary is a 
thermodynamical necessity in PZT system and further is 
very crucial for the high electromechanical coupling coef- 
ficients observed in MPB compositions. Two different ap- 
proaches for modelling the two-phase structures of PZT 
ceramics with compositions corresponding to the MPB 
were identified. The first is a method were space group 
symmetry was decreased so that both local and average 
crystal symmetries were modelled, whereas the second 
method was based on the highest space group symme- 
try compatible with the observed Bragg reflections and 
where local distortions were taken into account by se- 
lecting an appropriate line shape. Recently published 
Cm + Cc and Cc+ Cm structure models were compared 
with Cm + i?3c model where anisotropic line shapes were 
used. Whereas Cm + R3c model is a higher symmetry 
version of Cm + Cc model, Cc + Cm model was shown 



15 



to be inapplicable. Once the local distortions, resulting 
in anisotropic line broadening, are correctly taken into 
account, we do not find support for Cm + Cc model. In- 
stead we found that Cm + Ric model describes all Bragg 
peaks and their intensities well. It was further shown 
that the Cm + R3c model is consistent with the struc- 
tural features observed at other compositions and tem- 
peratures, which is particularly important in the vicinity 
of the phase boundary. Also issues related to the effects 
of a finite sample size used in transmission electron mi- 
croscopy and electron diffraction (ED) studies were dis- 
cussed and it was pointed out that while locally ordered 
areas may occur in nanosize samples, resulting in addi- 
tional reflections in ED patterns, they frequently do not 
appear in bulk samples. Thus, the same structural model 
might not be valid for nanoscale and bulk samples. It was 



further pointed out that these features are not limited to 
MPB compositions. Spatial variation of Pb ions shifts 
follows from the spatial variation in the i?-cation distri- 
bution. This is believed to be important for explaining 
the experimental structural data and for understanding 
electrical conductivity and piezoelectric properties. Con- 
nection between the structural parameters and piezoelec- 
tric properties were discussed in light of the recent com- 
putational studies. 
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